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Length-scale estimates for the LANS-a equations 
in terms of the Reynolds number 

J. D. Gibbon and D. D. Holm 
Department of Mathematics, Imperial College London, London SW7 2AZ, UK . 

Abstract 

Foias, Holm & Titi have settled the problem of existence and uniqueness for the 3D 
LANS-a equations on periodic box [0, L]^. There still remains the problem, first introduced 
by Doering and Foias [30] for the Navier-Stokes equations, of obtaining estimates in terms 
of the Reynolds number Re, whose character depends on the fluid response, as opposed to 
the Grashof number, whose character depends on the forcing. Re is defined as Re — Ul/v 
where [/ is a bounded spatio-temporally averaged Navier-Stokes velocity field and I the 
characteristic scale of the forcing. It is found that the inverse Kolmogorov length is estimated 
by (-X^^ < c{l/aY/^Re^/^. Moreover, the estimate of Foias, Holm & Titi for the fractal 
dimension of the global attractor, in terms of i?e, comes out to be 

where Va — (£/(i!Q;)^/^)^ and Vi — {L/tf'. It is also shown that there exists a series of time- 
averaged inverse squared length scales whose members, q), are estimated as {n > 1) 

^'^(k.I.o) < c„,aVcf^i?eT-ir(lni?e)" cii?e(ln i?e) . 

The upper bound on the first member of the hierarchy coincides with the inverse 

squared Taylor micro-scale to within log-corrections. 

1 Introduction 

1.1 Background to the LANS-a model 

Turbulence is a state of continual unrest, which arises as a fluid's response to forcing by stirring, 
for example, or by flowing along a wall. The scaling properties of turbulence are characterized 
by two dimensionless numbers. These are the Grashof number Gr, which measures forcing, and 
the Reynolds number Re which characterizes the fluid's response to the forcing. The turbulent 
response to forcing produces strong fluctuations in the fluid motion whose statistics obey power 
law spectra extending over a large range of length scales and time scales • This fluctuating 
multi-scale response is the hallmark of turbulence. Turbulence researchers often characterize 
the development of the multi-scale response as a "cascade" of kinetic energy rushing downward 
from the larger fluid motions due to forcing to the smaller and smaller circulations of eddies, 
sheets, and tubes of vorticity. In stretching themselves into finer and finer shapes, these vortical 
structures comprise the "sinews" of turbulence. 

Characteristic features of turbulence - its distribution of eddy sizes, shapes, speeds, vorticity, 
circulation and viscous dissipation - may all be captured by using the exact Navier-Stokes 
equations. These correctly predict how the cascade of turbulent kinetic energy and vorticity 
accelerates and how the sinews of turbulence stretch themselves into finer and finer scales, until 
their motions reach only a few molecular mean free paths, where they may finally be dissipated 
by viscosity into heat : for more details, see [H OlIS IHl El ■ However, the fidelity of the Navier- 
Stokes equations in capturing the cascade of turbulence is also the cause of serious problems in 
direct numerical simulations of turbulence. 

Based on Landau's heuristic idea that the number of active degrees of freedom required to 
simulate the turbulent cascade in high-Reynolds-number flows increases as it is clear 
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that this geometric rate of increase quickly outstrips the numerical resolution capabilities of 
even the largest computer. To make turbulence computable, scientists have developed various 
approximate models that halt the cascade into smaller, faster eddies. In most models, this effect 
is accomplished by causing the eddies below a certain size to dissipate computationally into 
heat. This dissipative imperative causes errors, however, because it damps out the variability 
(known as intermittency) in the larger-scale flow, which is caused by the myriad of small scales of 
motion interacting nonlinearly together in the fields of the larger motions. Thus, computational 
turbulence closure models based on enhancing viscous dissipation over its physical Navier-Stokes 
value run the risk of producing unrealistically low variability. 

Perhaps surprisingly, one of the first clues in understanding how to develop turbulence closure 
models without enhancing viscous dissipation came from the great mathematical analyst Leray 
[H] who showed how to regularize the Navier-Stokes equations by modifying their nonlinearity 
to the well-known form 

Vt + u-Vv + Vp = uAv + f{x) , diYV = 0, (1.1) 

with V = on the boundary. Here v is the (constant) kinematic viscosity coefficient, f{x) is the 
prescribed external force and u = Gg * v is a filtered version of the regularized velocity v. The 
filtering operation is defined hy Gs * v = J Gs{x, y)v{y) d^y for a symmetric kernel Gs{x, y) 
of characteristic width 5. The Navier-Stokes equations for v are recovered in the limit as 5 — > 0, 
so that u ^ V. For a review of the Leray regularization of the Navier-Stokes equations, see 
Gallavotti 

One of the points made in # is that the Leray regularization of the Navier-Stokes equations 
no longer satisfies the Kelvin circulation theorem, since for these equations 

^ (f V ■ dx = (f (vt + u-Vv + Vu^ • • dx / (f (lyAv + A ■ dx . (1.2) 
dt Jr{u) Jt{u) ^ ' Jt{u) ^ ' 

Remarkably, combining the process of Lagrangian averaging with Taylor's hypothesis (that 
fluctuations have such low power that they may be regarded as being carried along by the mean 
flow) leads to a regularized set of equations which answers the challenge of jH] to produce a 
regularization of the Navier-Stokes equations which do possess a Kelvin circulation theorem. 
These regularized equations comprise the LANS-a model^ 

Vt^u-Vv^SlviF -v^VP = vl\v^ j{x), divu = 0, (1.3) 

where P = p — iV(|ii|^ -|-a^|Vup). An equivalent alternative formulation is to rewrite ()1.3() as 

Vt — u X curl V = vAv — Vp + f{x) (1-4) 

where p = P + u- v. Usually, a is taken as a constant with dimension of length and the filtering 
relation u = Ga * v for the advection velocity in the LANS— a model is specified as 

v = u — a^Au. (1-5) 

Some remarks are in order here: 

1. The filtering kernel Ga for the LANS— a model thus turns out to be the Green's function 
for the Helmholtz operator (1 — a^A). 

^LANS is an acronym standing for "Lagrangian-averaged Navier-Stokes" while alpha (a) is the coherence 
length of the Lagrangian statistics; this term will be used hereafter. 
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2. As expected, the LANS— a motion equation satisfies the Kelvin circulation theorem: 

— (f V ■ dx = J) (vt + u ■ Vv + Vu^ ■ v) ■ dx = (f (i^Av + /) ■ dx . (1-6) 
dt Jr{u) Jr{u) ^ ^ Jr{u) ^ ' 

The circulation theorem tells us that the rate of change of circulation of momentum per 
unit mass v around a closed material loop r(w) moving with velocity u = Ga *v \s due to 
the integral around that loop of the tangential component of the sum over forces (viscous 
and external) acting on the fluid. 

3. This statement of the circulation theorem is also a mnemonic for deriving other regularized 
turbulence models of LANS— a type by specifying a different filtering kernel Ga jlOj - 

We have seen that the LANS— a model can be immediately derived from its circulation theorem. 
However, the approach used historically in Chen et al. 11 and Holm and Titi J2] for deriving the 
closed Eulerian form H1.3j) of the LANS— a motion equation was based on the combination of two 
other earlier results. First, the Lagrangian-averaged variational principle of Gjaja and Holm jl3j 
was applied for deriving the inviscid averaged nonlinear fluid equations, obtained by averaging 
Hamilton's principle for fluids over the rapid phase of their small turbulent circulations at flxed 
Lagrangian coordinate : this step had its own precedent in earlier work on Lagrangian-averaged 
fluid equations by Andrews and Mclntyre Second, the Euler-Poincare theory for contin- 

uum mechanics of Holm, Marsden and Ratiu was used for handling the Eulerian form of the 
resulting Lagrangian-averaged fluid variational principle. This second step determined the rela- 
tionships among the momentum per unit mass v, the Lagrangian-averaged velocity of the fluid 
u and the Lagrangian fluctuation statistics. Next, Taylor's hypothesis of frozen-in turbulence 
circulations was invoked for closing the Eulerian system of Lagrangian-averaged fluid equations, 
by obtaining the explicit relation v = u — a^Au. Finally, Navier-Stokes Eulerian viscous dissipa- 
tion was added, so that viscosity would cause diffusion of the newly-deflned Lagrangian-average 
momentum and monotonic decrease of its total Lagrangian-averaged energy. 

At this point, one may regard the LANS— q model as an alternative regularization of the 
Navier-Stokes equations and re-examine its properties from the viewpoint of Leray's analysis: in 
fact the problem of existence and uniqueness for the LANS-a equations has already been settled 
by Foias, Holm & Titi FT6l. In addition, the same ideas which restore Kelvin's circulation 
theorem to Leray's regularization of the Navier-Stokes equations also turn out to provide a basis 
for deriving candidate equations for Large Eddy Simulations (LES) of turbulence. Conversely, 
other proposed LES models of turbulence may lead to likely candidates for application of Leray's 
analysis. In this way, the classical Leray analysis of the Navier-Stokes equations flnds itself a 
new role in the study of the analytical properties of turbulence models. Indeed, the Leray model 
itself was recently found to be a viable candidate for LES computational modeling of turbulence 

[naini. 

1.2 Estimates in terms of Reynolds number 

Motivated by the Navier-Stokes equations, in an early and progressive paper Ruelle J^, discussed 
how ideas in dynamical systems might be extended to the inflnite dimensional case by counting 
the number of positive characteristic exponents. Certain general flniteness assumptions were 
made about the nature of the Navier-Stokes equations without formally using the concept of 
a global attractor. These ideas were taken to another level by Constantin and Foias [T^ who 
used the idea of a global attractor A for a partial differential equation to determine the number 
of growing Lyapunov exponents and thence, through a rigorous generalization of the Kaplan- 
Yorke formula, they were able to get formal upper bounds on dL{A), the Lyapunov dimension 
of A, which itself bounds above dpiA) and dniA). Importantly these bounds depend only upon 
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time-averages and thus only on the long-time dynamics on A; see also ^20^. When applied to the 
two-dimensional Navier-Stokes equations, for which a global attractor exists, these methods gave 
good, sharp estimates to within logarithmic corrections [25 when used in conjunction with an 
L°° estimate of Constantin j^. The development of the key features of these ideas can be found 
in |2l EHI , including the role played by the inequalities of Lieb and Thirring |23] . Apart from 
the two-dimensional Navier-Stokes equations, several important partial differential equations 
possess a global attractor, such as the the two-dimensional complex Ginzburg-Landau equation 
and the one-dimensional Kuramoto-Sivashinky equation. Unfortunately, the three-dimensional 
Navier-Stokes equations are not among them; the key element is the lack of a proof of existence 
and uniqueness without which the existence of a global attractor remains open. 

In settling the question of existence and uniqueness for the LANS-a equations, Foias, Holm 
&: Titi were able to use the machinery developed for global attractors to find estimates for 
the dimension of its global attractor A (and other important quantities) in terms of the Grashof 
number Gr, which is a dimensionless control parameter dependent only on the ratio of the 
forcing to the viscosity u. Therefore, as a regularization of the three-dimensional Navier-Stokes 
equations with many similar features to its parent, the LANS-q equations equations possess this 
extra key property. What remains to be proved, however, is whether the estimates in \16^ can 
he evaluated in terms of the Reynolds number, whose character depends on the fluid response to 
the forcing, and which is intrinsically a property of Navier-Stokes solutions. The advantage of 
this further step lies in the fact that the engineering and physics communities express their ideas 
about turbulence in terms of the Reynolds number. It also allows us to make direct comparisons 
between estimates for the two equations where they exist (see Table 1). 

For simplicity the LANS-a equations ^\.\^ . or their alternative form in (jl.4|) . will be consid- 
ered on a periodic domain [0, L]^^^ with forcing f{x) taken to be L^-bounded of narrow-band 
type^ with a single length-scale i. (see [2131^132]) 

l|v"/||2«r"!|/||2. (1.7) 

With frms = -^"'^^^ll/lb, where and H/lH = \ f\^dV, the standard definition of the Grashof 
number in d-dimensions is 

Gr = llnni, (1.8) 

Analytical estimates in Navier-Stokes theory have traditionally been expressed in terms of Gr 
[3 IS EEl HZl I2HI IZH] but these are difficult to compare with the results of Kolmogorov scaling 
theories which are expressed in terms of Reynolds number 2 . A good definition of this is 

Ee = ^ C/2 = L-^(||«||2) (1.9) 

where (•) is the long-time-average 



1 



(5(-)) = limt^oo- / g{T)dT. (1.10) 







Doering and Foias j3L)j addressed this problem recently and have shown that in the limit Gr 
solutions of the (i-dimensional Navier-Stokes equations 



00, 



lit + w • Vw + Vf> = z^Aii + /(a;) , divt> = 0. (1.11) 

must satisfy^ 

Gr < c{Re^ ^ Re) . (1.12) 



^The restriction to narrow-band forcing can be relaxed at the cost of more parameters in the problem. 
■^This result is not advertised in |30| but it follows immediately from their equation (48). |3U| also contains 
another result that the energy dissipation rate e has a lower bound proportional to Gr: see Appendix 1A.2II . 
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Time-average 


NS 


LANS-a 


Eqn number 








(ITTT^I) 


(^X^^ (improved) 


unknown 


{£/ay/'^Re^/^ 


(ITHI) 


a^^i/-^ {H2) 


unknown 


VeRe^ 


(ESI) 




unknown 


F„y//'(L2Ai)-9/8 i?e9/4 






unknown 


VaRe^" {In Re)^ 






i?e(ln Re) 


Re{lnRe) 






unknown 


Re^y^ 




aiiy-^ (||Vw||oo) 


unknown 




(EH) 




y^i?e3-5/2"(lni?e)^ 


Va^"^ Re^"/^{lnRe)^ 





Table 1: Comparison of various time-average bounds for the Navier-Stokes and LANS-a equations with 
constants omitted. Note that in the cases where there is no known equivalent upper bound for the 
Navier-Stokes equations, the LANS-a upper bounds blow up as a 0: A„ is defined in (|1.23(l . Lines 5-8 
are a summary of the results of Theorem|21in 21 The Navier-Stokes estimate for £ (^,1,0) on the last line 
is there for comparative purposes and is taken from I82| . 

Using the relation in (|1.12|) , Doering and Gibbon |3H I32j have re-expressed some 3-D Navier- 
Stokes estimates in terms of Re (see Table 1). The problem, however, is less simple than 
substituting (|1.12|1 into standard results, although this works well enough for point-wise esti- 
mates jl6| . Time averages are more subtle and exploit how the average velocity U within Re is 
related to (^HiiHl). As an illustration, let us consider the Navier-Stokes equations (|1.1H) whose 
energy dissipation rate is e = v (||Vti||2) L~'^. Standard estimates show that its upper bound is 
proportional to Gr^. By (|1.12|) . this turns into Re^, which is not sharp. Now we estimate this 
a different way consider Leray's energy inequality 

hj^lMl < -^l|Vw||i + ||/||2||«||2. (1.13) 

Time-averaging p.l3() and using (|1.9() and ()1.12|) yields 

e < v^r^GrRe < cv^l'^ [Re^ + Re) , (1.14) 

which is a considerable improvement. To leading order the inverse Kolmogorov length A^^ = 
{e/v'^y^^ is then bounded above by 

£X-^<cRe^/^. (1.15) 

This estimate now conforms with the generally accepted scaling law for the inverse Kolmogorov 
length with the Reynolds number ^H]. 

The relation in (|1.12j) is essentially a Navier-Stokes result and thus needs re-proving for 
the LANS-a equations. In turns out to be true but the proof is a non-trivial extension of the 
method in [S^; the whole of Appendix (|A.1|) is devoted to this proof. As will be shown in fj^l 
the estimate for (Hi) in (|1.14j) can be improved to 

//-X 1/4 

e<ci^H-^a-^Re^/\ iX^^ < c i^j Re>'\ (1.16) 

Moreover, in ^the estimate by Foias, Holm and Titi ^H] (see also for the fractal dimension 
dF{.A) of the global attractor A is considered in the light of these i?e-bounds. Using their 
estimate in terms of the generalized a-dependent dissipation rate, we show that 
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where the two dimensionless volumes Vi and Va are defined by 

and Ai > is smallest eigenvalue of the Stokes operator. Given our definition of Re in ()1.9() . 
this iie^/^ estimate is consistent with scaling theories of turbulence but it does not survive in 
the Navier-Stokes limit because the volume Va blows up as a ^ 0. 

The Re^^'^ estimate for dp also gives an idea of how many degrees of freedom, in Landau's 
sense, exist in a turbulent flow - indeed this is exactly result predicted by Landau for the 
Navier-Stokes equations^. What is not taken into account in this picture is the effect of strong 
dissipation-range intermittency where significant energy lies in wave-numbers larger than A^^. 
In this case estimates are needed for length-scales that are associated with higher derivatives. 
This idea has been investigated in |3H I32j where estimates were found for time-averaged Navier- 
Stokes quantities. Those for the LANS-a equations should be much better because of their 
enhanced regularity properties |16j . In fQ] we combine the forcing with higher derivatives of the 
velocity field in the form 

F„ = i/„ + r2||V"/||i, (L19) 

where {Gr\nGr) ' is a characteristic time: see Appendix (|A.2|) . We also form the 

combination 

Jn = Fn + a^Fn+i (L20) 
and use it to define a set of inverse length scales, or time-dependent wave-numbers, 

K„,oW = (^)" • ^^-^^^ 

In the a — s- limit, the K^"g behave as the 2nth-moments of the energy spectrum. Theorem |21 
of ^proves that the time average of their squares must obey 

{nlfl) < Cn,aVa " iJe'^" (In i?e) " + ciRelnRe, (L22) 
for the LANS-a model, where 



Note that the n = 1 estimate in (|1.22|) scales with Re the same as the Taylor micro-scale. This 
best that could be achieved for the full 31? Navier-Stokes equations was [311 132j 

^(Kn,o) <CnVi>Re^'^-{\nRe)^ +ciRe\nRe, (1.24) 

which appears in the last line of Table 1. With the exponent of unity in the time-average -if 
indeed a solution exists at all - this is not only much worse that H1.22() but represents only weak 
solutions, as opposed to the strong solutions of Foias, Holm and Titi jl6j . The fact that no 
upper bound is known to exist for (k^q) Navier-Stokes equations is consistent with 

the fact that the dimensionless volume Va blows up as a —> 0. The k„^o could even become 
singular in this limit. 



^However, given the improved Re^^^ inverse Kolmogorov estimate in 11.1611 for the LANS-a model, it is possible 
that the sharp estimate for dp is proportional to Re^^^^ . 
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2 Properties of the product j^u ■ v dV 

Foias, Holm and Titi 16 noted that the product u ■ v has two convenient properties 

/ u-vdV= [ {\u\^ + a'^\Vu\^}dV (2.1) 
Jn Jn 



and 



— I u ■ V dV 

dt 



/ {ut ■ V + u ■ vt) dV 
Jn 

= {ut ■ {I - cPA)u + u ■ Vt)dV 
Jn 

= / {u ■ [l - a^A)ut] +u ■ vt} dV 
Jn 

= 2 u-vtdV (2.2) 
Jn 

where two integrations by parts have occurred between the second and third lines. Now define 

Hn= [ |V"-updy= [ \curru\'^dV (2.3) 
Jn Jn 

this being true on a periodic domain because divu = 0. Clearly we have the bound 

7 p 

\—(Ho + a^Hi) = -u[Hi+Q^H2)+ u-fdV 
dt Jn 

,2 



< -iy{Hi + a'H2) + \\u\\2\\fh. (2.4) 



One can then calculate an absorbing ball for Hi with ease (see |16j). It is also possible to 
estimate the time averages {Hi) and {H2) which can be found in the same manner as in (|1.13|) 
to satisfy 

uL-'^ {Hi + a^i^s) < v^r^Re Gr < c v^T^Re^ . (2.5) 
The upper bound on (-ff2), written as 

a^iv-'^ {H2) <cViRe^ (2.6) 

can then be used to improve the estimate for {Hi) by using the fact that {Hi) < {Hq)^^'^ (^2)"'^^^ 
and that U'^ = L'^ {Hq). We find that 

{Hi) < cu^L^r^a-^Re^/^ , (2.7) 

and so 



a 



Hence the energy dissipation rate e is also bounded above by Re^l"^ but the improved estimate 
blows up when a — > 0; no equivalent result is implied for the 31? Navier-Stokes equations. 

3 A i?e^/^ bound for the attractor dimension 

Foias, Holm and Titi ^B] made two independent estimates of the fractal dimension dp{A) of the 
global attractor A. The first was in terms of Gr but the second estimate was made in terms of 
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the "energy dissipation rate" e; this phrase has been put in inverted commas because it includes 
the iif2-norm, whereas conventionally only the i^i-norm is used. Their definition of e is 

€ = Xl^\{Hi + a^H2) (3.1) 

where Ai is the smallest eigenvalue of the Stokes' operator which has the dimension of an inverse 
length squared. Foias, Holm and Titi then proved that 

> -3/2 / _ X 3/4 

dF{A)<c^^^(^A . (3.2) 



(a2Ai)3/4 

They then defined a Kolmogorov length as = (e/i/'^)^^^ which turns (|3.2|) into 

dpiA) < C—r^—rj- • (3.3) 

(a2Ai)3/4 (4a}/^)3 

The problem here lies in interpretation: is not the conventional Kolmogorov length because e 
is not the Navier-Stokes energy dissipation rate e = {Hi) L~^. Instead we take an alternative 
route and use the estimate for i^Hi + H2) from (|2.5() . which we repeat here 

uL-^ {Hi + a^H2) < c u^r^Re^ . (3.4) 

Thus 

e<c(LXyyu^r^Re\ (3.5) 
which turns the result of Foias, Holm and Titi |16j into 

V v^l'^ 

^^^-^^^^(iS^^^'^'' ^'-'^ 

where l?\i = Att'^. The right hand side blows up as a ^ through V^. Despite this, the Re^/^ 
estimate is, to our belief, the first time this has been achieved with this definition of Re. As 
has often been pointed out, this upper bound is also valid for the Hausdorff dimension dniA) 
because dniA) < dpiA). 



4 A Theorem involving higher derivatives 

In terms of the number of degrees of freedom, the result in H3.6() says that i?e3/4 X i?e3/4 x i^e^/^ 
resolution grid points are needed. However, as explained in ^1.21 what is not taken into account 
in attractor dimension estimates is the effect of strong dissipation-range intermittency where 
significant energy lies in wave-numbers larger than A^^. In this case estimates are needed for 
length-scales that are associated with higher derivatives. To obtain such estimates we begin by 
forming the combination 

F„ = i/„ + r2||V"/||i, (4.1) 

where the quantity r 

T = fu-^{Gr In Gr)-^/^ (4.2) 

where the InGr-term is there for reasons explained in Appendix ()A.2|) . We also define the 
combination 

Jn=Fn + a^Fn+i . (4.3) 

The ultimate aim is to find time-averaged estimates for the Kn^r that appeared in 1)1.21(1 : this 
will be the subject of the next section. In preparation, we prove the following result: 
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Theorem 1 As Gr oo, for n > 1, 1 < p < n, Jn satisfies 



i+i 

< — ^ It 

dt 



J. 



n—p 



and, for n = 0, 



1 dJp 
dt 



(4.4) 



(4.5) 



Proof: The results on the pairing of Hq and Hi in l|2.1j) and (|2.2j) apply more generally: 

/" (V"u) • (V%) = [ {\V''u\^ + a^\V'^+^u\^}dV = H,, + a'^Hn+i. (4.6) 
Jn Jn 



Thus 



^^(Hn + a^Hn+i) = [ (V"ii) • (V"^t) dV 
dt 

which, from l\lA\i . gives the estimate 



[ (V"m) • {z^AV"-y - V"(-u X curlt;)} 



< -U [Hn+l + + 



[ (V"+^u) (V"-^^* X curli;)) dV 
Jn 



(4.7) 



(4.8) 



(4.9) 



Upon separating the two constituent parts oi v = u — a'^Au within the last term in H4.8p . the 
first is found to satisfy 



/ (V"+^u) • (V"-^(ix X curl-u)) dV 
Jn 



^ II "||oo-nn+l-"n 



E C^-'W'^'^uUV^uW,, (4.10) 



n-1 



where p, q must satisfy p ^ + q ^ = 1/2 according to Holder's inequality. The first term on the 
RHS of (|4.1flj) is the m = term. Now we use the two Gagliardo-Nirenberg inequalities 



where 



iV^ullp < cIlV^-ull^^llull^"" 



1 m /In 
- - ^ = am (^2 " d 



iV^-^-ujlg < c||V"'u||^■'"||-u||^^" 



l n — m 



q d 



1 n 
2~ d 



(4.11) 



(4.12) 



Adding the two and noting that p + q =5 implies that am + = 1- In fact this is true in 
d-dimensions. Thus the last term in (|4.1fl|) is estimated by 



n-1 



(1) f7l/2H-l/2 



(4.13) 



in=l 



We may approach the second constituent part of v in the same manner as 1)4. lOp . After an 
integration by parts we find 



/ (V"+2ti) . (V"-2 [u X curl(-a2Aw)) dV 
Jn 



^ (2) 2r7l/2 rrl/2 ii i 



(4.14) 
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The estimate ()4.8p now beomes 
d 



+ \\V-fhHll' + a'\\V-+'fhHT 



+1 • 



(4.15) 



To turn this into an inequahty for the we add and subtract to the negative terms on the 
right hand side, and break up the last terms to form the group of terms designated as 



Xn = WV^fhHl'^ + a^ll V^+Vlb^^fi + |i^rV||V"/||2 + ii.rVa2||v"+i/i (4.16) 



with 



d 
^dt 
where 



1- (F„ + a^Fn+i) < -\v {Fn+i + a^Fn+2) + Cn (F„ + Fn+i) + Xn (4.17) 



Xn ^ S —H, 



2gT^ 2P J 



r1|V"/||2K<!|^n+l + 



+ — I r^||V"+VI 



To make the coefficients of Hn and t ||V"/||2 equal, choose g to satisfy 

9 -^-72=0- 



That is 



2^ = ^ + 



\2 4 
£2j 



1/2 



Given r in 1)4.2(1 we have g ^ t ^ as Gr — > 00. Consequently, 

r"^ = fu-^{GrlnGr)^2 < c f Re{ln Re) 

and 

Xn < (F„ + a^Fn+i) . 
When applied to (|4.17() . the previous result yields 

which is 

The following Lemma deals with the — Jn+i-term in ()4.23|) : 
Lemma 1 For 1 < p < n, the Jn satisfy 

1 T 

Jn+l ^ 2 jl/p ■ 

Proof: Firstly from ()4.1|) . by writing Fn in Fourier transforms one finds 



Fn = Hn + T^V''f\\l 



(4.18) 

(4.19) 
(4.20) 

(4.21) 
(4.22) 

(4.23) 



(4.24) 



(4.25) 



where includes the three components of the fluid velocity and the three components of the 
forcing. Using Holder's inequality produces 



^2(n+g)[^^|2y+. ( ^2(n-p)j^|2^ ^y^ < _ 



(4.26) 
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Another application of Holder's inequality also gives the standard result 

/TV \ P+1 N 



(4.27) 



,1=1 / 1=1 



Now we find upper and lower bounds for the combination Fi^^ + a'^^^'^F^^l. Inequality (|4.26|) 
is used to find an upper bound 

-r„ -r UL ^n+1 — ^n+l^n~p ^ « ^n+2^n+l-p 

< {Fn+1 + a^Fn+2Y {Fn-p + Fn+l-p) . (4.28) 

A lower bound comes from inequality (|4.27|) with N = 2 

FP+' + [a'Fn+ir+' > {Fn + a'Fn+iY^' . (4.29) 

Thus we have 

< Jl^,Jn-p (4.30) 
which gives the result. □ 
The n = result 1)4. 5p follows from 1)2. 4p by using the same methods. ■ 

5 Length scales 

Now define the quantities r < n (we take r = n — p) 

1 

(T \ 2(n-r) 

which act as t-dependent wave-numbers and thus have the dimension of inverse length scales. 
For r = they are analogous to the 2n-th moments of the energy spectrum. In what follows we 
find upper bounds for (k^ q). These are bounds on an infinite series of inverse squared length 
scales. 

In Theorem 1 for r < n, dividing by the J„ & time averaging , we have 

(k?,o> < cii-'^RelnRe, (5.2) 

< Cn,aiy'^ {Mil,) +cir^Re\nRe. (5.3) 

The results for this can be summarized in the following: 
Theorem 2 As Gr oo, the four time averages are estimated as 

n — l 11 7 1 

< Cn,aVa'' Re^-^^ {In Re)-^ +ciReln Re, (5.4) 

^"^{i^io) < ciRelnRe, (5.5) 

i'^^-'iMlo) < cV.Re''/\ (5.6) 

aii^~' {\\Vu\U < cyi/V//'i?e35A6. (5.7) 

Remark: These four results are also listed in the Table of ^1.21 
Proof: Agmon's inequality says that 

||w||L < c ||V2-u||2||Vu||2 = cHI''^h\''^ . (5.8) 
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Moreover, we also notice that because Hi < H^^'^ H^"^ and U = L {Hq)^^"^ we have 

{Mlf < c {Ho)'/' {H.f < cL^'^U {H^f . (5.9) 

which simphfies to^ 

f^-'{\\u\\l)<cVaRe''/^ (5.10) 
which is H5.6|) . With the dimensionless volumes Va and Vi defined in p.l8|l . we have 

^^(k?_o> < ciRelnRe, (5.11) 

i^i^l^r) < Cn,aVaRe''/'' + CiRelnRe. (5.12) 

For r = we can improve this by writing 



(('^n,l)""'('^l,o)") 



< «i> " ('^lo)" , (5-13) 
and then using the estimates in (|5.11j) and (|5.12|) for n > 1, which gives 

n — 1 11 7 1 

(k^,o) < Cn.aVa " Re~ (in Re)"^ + CiRelnRe. (5.14) 

This is (j5.4j) . Note that when n = 1 we return to < ciRelnRe. 

It is also possible to estimate (j|Vii||oo) 

(llVtzlU) < c {hI/'hI^') = c {.li^H',/') < c {4,,f" ■ (5.15) 

Using ((231) and (I^TT^ we find 

alu-^ (llV^lloo) < c [hI'^hI'^) < cyi/V//'i?e35/i6 , (5.16) 

which is ((^ . ■ 
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nating discussions of the LANS-a equations and their analytical properties. We also thank C. 
R. Doering for encouraging explanations of the results of 30 . The work of DDH was partially 
supported by US DOE, Office of Science, Applied Mathematics program of the Mathematical, 
Information, and Computational Sciences Division (MICS). 

A Issues concerning the forcing 
A.l Bounds concerning Gr and Re 

Doering and Foias [30] split the forcing function f{x) into its magnitude F and its "shape" $ 
such that 

f{x) = F^{r'x) (A.l) 

^Herein lies the difference between 3-D LANS-a equations and the 3-D Navier-Stokes equations: in the former 
we have a bound on {-ff2) from 12.51 which is missing in the latter. 
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where i is the longest length scale in the force. On the unit torus I^, $ is a mean-zero, divergence- 
free vector field with the chosen normalization property 

/ |V-^*|'d'^y = l. (A.2) 

L^-norms of / on are 

||V^/|1^ = CNr^^L'^F^ (A.3) 
where the coefficients C n refer to the shape of the force but not its magnitude 

CM = 5]|27rn|2^|*„|2. (A.4) 

n 

Doering and Foias j3Uj showed that various bounds exist such as (among others) 

||VA-^Vlloc = I)Afi^^''''-^ (A.5) 

The energy dissipation rate e is 

e = (^i^L-'^ j I Vup dV^ = uL-'^ {Hi) . (A.6) 

In terms of F the Grashof number in (jl.Tf) becomes 

Gr = Fl^jv^ (A.7) 



and the Taylor micro-scale At is related to IJ via At = ^/vV^ /e , which is consistent with the 
definition A^^ = {H/) / {Hq). 

The LANS-a equations ()1.3|) and ()1.4|) can also formally be re-written as 

Ut + u Vu- uAu + Vp = (1 - a^A)"^ {f{x) + a^divT} (A.8) 
where the tensor T is defined as 

T = Vii • Vit + Vii • Vii^ - Vii^ • Vu . (A.9) 
Following the procedure in |^ (pg 296 equation (2.9)) and multiplying by (— A~*^)/ we have 

/ u-[{-A~'')f]dV = -V I Au.[{-A-^')f]- f u-Vu-[{-A-'')f]dV 

+ [ [i-A'^')f] ■ (1 - q2a)-i {/ + Q^divT} dV . (A.IO) 

Now there are two strategies: 

1) To prove that Gr < cRe^: integrate all the terms by parts, and take the time average 



d_ 
Jt 



L-'' I V~^'f ■{l-a'Ay^V-^'fdv) < I u-[{-A-^"+^)f\dV 

Jid / \ Jid I 

- Il-^ ^^u-m-A-'^'m-udv^ (A.ii) 



+ Ih-'^l |V[(-A-^'^)/]| 



:i -a^A) 



-T 



dV 
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Using the scaling properties of $ and a Fourier transform on the last term 



L-'^ [ |V[(-A-^^)/]| 



a 



-T 



21,2 



1 + a^P 



\u\^dVk 



where ()A.5|) defines Dm- Thus ()A.11|) turns into 



n2 «2M 



Co 



1 + a 



20-2 



<cii^Ff''-'U + C2f^-'FU\ 



(A.12) 



(A.13) 



where the C/^-term contains the contributions from both nonlinear terms and the constants (not 
explicitly given) contain the shape of the body forcing. Using (|A.7|) . HA.13|) becomes 



Gr<c {Re + Re^) 



Gr — > oo . 



(A.14) 



the only difference from the Navier-Stokes equations being the value of the constant. 

2) To prove that e > cv^l^^L^^Gr: return to HA.10|) and take a different route. Firstly in 
the Laplacian term use one derivative on u and another on the forcing. Then keep the u ■ Vu 
advection term. Finally integrate the T-term by parts and exploit the fact that (1 — a^A)~^ is 
a symmetric operator. 



+ U 



+ {L 



u\ |Vit| 

a" 



)f\dV 



(A.15) 



V[(-A-^-^)/] 



(l-a2A) 



ITI dV 



Thus (fIT5|l turns into 



F 



2n2M 



Co 



1 + a 



20-2 



< CiFz.l/2^2M-l^l/2 ^ ^^^^-l/2^2Af^l/2^ ^ Csa^l^^'-^ Fl^'h . 



(A.16) 



We assume the ordering of length scales as L > £ > a > Xt and use the fact that U"^ = Xl 
L^eu~^. Then HA.16|) becomes 

|coF < ciu^/^rh^/^ + (caL + c^a^r^) < ciu^^^rh^/^ + c^i^-^Le . (A.17) 

For Gr — > oo the last term on the RHS is the dominant one: we have 

e > cs^^^L^V^Gr (A.18) 

which, with a different constant, agrees with the result in Doering and Foias [HHj. This inequality 
is used in the next subsection. 



A. 2 Forcing & the fluid response 

For technical reasons, we must address the possibility that in their evolution the quantities Hn 
might take small values. Thus we need to circumvent problems that may arise when dividing by 
these (squared) semi-norms. We follow Doering and Gibbon j31j who introduced the modified 
quantities in 1)4. 1() 

Fn = Hn + T^V^f\\l (A.19) 
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where the "time-scale" r is to be chosen for our convenience. So long as r 7^ 0, the -F„ are 
bounded away from zero by the explicit value r^L^^"^*^/^^^. Moreover, we may choose r to 
depend on the parameters of the problem such that ~ {Hn) as Gr 00. To see how to 
achieve this, let us define 

T = fu-^ {Gr In Gr)-^/^ . (A.20) 
Then the additional term in ()A.19|) is 

r^llV^/lli = L'u-h'-'y?^,{Gr\nGr)-' 

= u'^r^'^''+^^L^Gr{\nGr)-^ (A.21) 

Recalling the a priori bound on the far right hand side of HA.18|) 

(2n-l) 



C6 



(^j) L-2(-i)(/?i)(lnGr)-i (A.22) 



Using Poincare's inequality in the form Hi < (27rL)^(" ^^Hn, as Gr ^ 00 we have 

2||v7nf||2 /j^\{2n-l) 

^M<,.(^) (1„G,.)-' (A.23) 

Hence, the additional forcing term in ()A.19() becomes negligible with respect to (Hn) as Gr 00, 
so the forcing does not dominate the response. 
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